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ABSTRACT 

Using the clathrate hydrates trapping theory, we discuss the enrichments 
in volatiles in the atmosphere of Jupiter measured by the Galileo probe in the 
framework of new extended core-accretion planet formation models including mi- 
gration and disk evolution. We construct a self-consistent model in which the 
volatile content of planetesimals accreted during the formation of Jupiter is calcu- 
lated from the thermodynamical evolution of the disk. Assuming C02:CO:CH4 = 
30:10:1 (ratios compatible with ISM measurements), we show that we can explain 
the enrichments in volatiles in a way compatible with the recent constraints set 
from internal structure modeling on the total amount of heavy elements present 
in the planet. 



1. Introduction 

The formation and structure of giant planets, in particular those of the solar system for 
which numerous observational constraints exist, has been the subject of many studies during 
the recent years. In the present paper, we show that the core accretion scenario (Pollack et 
al. 1996, hereafter referred to as P96) extended to include migration and disk evolution can 
actually lead to the formation of a Jupiter-like planet which accounts for both the constraints 
set by internal structure models of this planet and the measurements of the abundances of 
volatiles in its atmosphere. 

Detailed internal structure models of Jupiter matching many observational constraints 
(mass, radius, surface temperature, gravitational moments, etc.) have been constructed 
recently (Guillot at al. 2004, hereafter referred to as G04). From this modeling, the mass of 
the core of the planet and the total amount of heavy elements present can be constrained. 
G04 found that this maximum amount of heavy elements present in Jupiter is of the order of 
~ 42M e , whereas the mass of the core can vary anywhere from to 13M®, the uncertainties 
being essentially due to our poor knowledge of the equation of state of hydrogen. Note that 
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the inferred present day core mass can differ significantly from the one at the end of the 
formation process depending upon the extent of core erosion (G04). The total amount of 
heavy elements present, on the other hand, is not affected by this process. 

Abundances of volatile species in Jupiter's atmosphere have been measured using the 
mass spectrometer on-board the Galileo probe (Atreya et al. 1999, Mahaffy et al. 2000). 
These measurements reveal that the giant planet's atmosphere is enriched by a factor of 
~ 3 in Ar, Kr, Xe, C, N and S compared to the solar abundances of Anders & Grevesse 
(1989, hereafter AG89). In order to explain these enrichments, Owen et al. (1999) suggested 
that these volatile species were trapped in planetesimals made of amorphous ice which were 
accreted by Jupiter. Given the low temperature required for the existence of amorphous 
ice, this suggestion implies one of the following three possibilities: 1) Jupiter was formed 
in sufficiently cold regions (about 30 AU) and subsequently migrated inwards to its present 
location, or 2) the solar nebula was substantially cooler at 5 AU than currently imagined, 
or 3) Jupiter accreted a substantial amount of ices originating from distant regions. 

Another explanation has been proposed by Gautier et al. (2001a,b) and Hersant et al 
(2004) who suggested that volatile species can also be trapped by H 2 ice at much higher 
temperature in forms of clathrate hydrates or hydrates. Hence, the measured enrichments 
could also be explained by the accretion of icy planetesimals originating from the local, 
warmer region of the nebula. However, the large amount of water ice required to trap the 
volatile species in these models implies a total heavy element content in Jupiter incompatible 
with the constraints derived by G04. Finally, while these models are based on the P96 
calculations, they are not fully consistent with these calculations for a number of reasons. 
First, they assume that planetesimals are captured by Jupiter during the late stages of the 
formation process, while in the models of P96, ~ 10M e of heavy elements are accreted during 
the early phases. Second, the thermodynamical conditions of the disk evolution models used 
in Gautier et al. (2001a,b) and Hersant et al. (2004) differ from the ones used as outer 
boundary conditions for the planet formation calculations of P96. 

The standard giant planets formation model by P96 has been recently extended by Al- 
ibert et al. (2004,2005a) who included migration of the growing planet and disk evolution. 
The authors show that these two phenomena can have significant consequences on the for- 
mation process. This motivated us to recalculate the enrichments in volatiles in Jupiter in a 
way fully consistent with these new formation models. For this, we use the calculations by 
Alibert et al. (2005b, hereafter referred to as A05) to compute the thermodynamical prop- 
erties of the disk, which allow the composition of the ices to be determined, and the time 
dependent accretion rate of planetesimal during the growth of the planet. As we will show, 
this fully self-consistent approach can explain the measured enrichments while requiring a 
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total abundance of heavy elements in agreement with the constraints derived by G04. 

2. Theoretical models 

2.1. Formation of Jupiter 

The model of the formation of Jupiter we consider in this work is the one presented in 
A05 which is based on the extended core-accretion formation approach described by Alibert 
et al. (2004,2005a). In this model, the giant planet forms from an embryo originally located 
at ~ 9 AU 1 which migrates inwards and stops at the current position of Jupiter at the time 
the disk disappears. 

This model takes into account the evolution of the protoplanetary disk in the framework 
of the a formalism (Shakura & Sunyaev 1973). The thermodynamical properties of the disk 
as function of position and time (T, P and E) are calculated by solving the vertical structure 
equations using the method presented in Alibert et al. (2005a). These quantities are used to 
determine the composition of the ices incorporated in the planetesimals (see next Section). 

2.2. Composition of planetesimals in Jupiter's feeding zone 

We assume that volatiles have been trapped during the cooling of the nebula in plan- 
etesimals either in form of pure condensates or, following Gautier et al. (2001a,b), in form 
of hydrates or clathrate hydrates. Once condensed, ices are assumed to decouple from the 
gas, to be incorporated into growing planetesimals which may subsequently be accreted by 
the forming Jupiter. 

Figure 1 shows cooling curves of the nebula at 5 and 15 AU derived from the A05 model, 
as well as the condensation curves for the various ices considered in this work. The stability 
curves of clathrate hydrates are derived from Lunine & Stevenson (1985) whereas for the 
pure C0 2 condensate a fit to the existing experimental data (Lide 1999, pp 6-59) is used. 
This figure provides the condensation sequence of the different volatiles initially existing in 
vapor phase inside Jupiter's feeding zone and the thermodynamical conditions (T, P and S) 
at which the different ices are formed correspond to the intersection between the cooling 
curve and the stability curve of the different condensates. 



1 The precise starting location of the embryo depends upon the assumed migration rate and can range 
from 9 to 14 AU, see A05. 
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In this work, we assume that the amount of water available along the migration pathway 
of the forming Jupiter is at least large enough to trap all the volatile species. With our 
assumption of decoupling upon condensation, this implies a minimum abundance of water 
relative to H 2 necessary of 

clement i H^,^H 2 0,-lH 2 0) 

where x« is the molar mixing ratio of the volatile % in vapor phase with respect to H 2 in the 
solar nebula, 7$ is the required number of water molecules to form the corresponding hydrate 
or clathrate hydrate (see Lunine & Stevenson 1985; 7$ = in the case of pure condensate). 
S(i?, P u Tj) is the surface density of the solar nebula at the distance R from the Sun, at the 
intersection of the cooling curve and the condensation curve of the species i. 

Finally, for a given abundance of water, we estimate Yj, the volatile % to water mass 
ratio in icy planetesimals, by the relation given by Mousis & Gautier (2004): 

Y _ = X i T,(R,Pi,Tj) ^ 
Xh 2 o S(i?, Ph 2 o, Th 2 o) 

where Xi and Xh 2 o are respectively the mass mixing ratio of the volatile i and of H 2 with 
respect to H 2 , in the vapor phase of the solar nebula. 



3. Results 

3.1. Initial ratios of C0 2 :CO:CH 4 and N 2 :NH 3 

In order to be able to compute the composition of the planetesimals, a number of 
assumptions have to be made. In the present paper, we assume that in the solar nebula gas 
phase the elements are in solar abundance (AG89, see Table 1). Note that we also assume 
that the amount of water ice is sufficient to trap all volatiles. As we shall see in Section 
3.2, this corresponds to an abundance of water relative to H 2 in vapor phase in excess of 
the value derived from solar abundances (see Table 1). Sedimentation and drift due to gas 
drag on the icy grains could explain this overabundance of water (Supulver and Lin 2000). 
Furthermore, we assume that 0,C and N exist only under the form of H 2 0, C0 2 , CO, CH 4 , 
N 2 and NH 3 . S is present only in the form of H 2 S and other sulfur components (see next 
Section). 

Starting with ISM initial values, Mousis et al. (2002) have calculated the evolution of 
the CO:CH 4 and N 2 :NH 3 ratios in vapor phase in the solar nebula by taking into account 
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turbulent diffusion and kinetics of chemical conversions. They demonstrated that these ratios 
remain quasi-identical to the ISM values. Therefore, we adopt CO:CH 4 = 10 in the vapor 
phase of the solar nebula, a value consistent with the ISM measurements of Allamandola et 
al. (1999). Moreover, observations of ISM ices reveal the presence of large amount of CO2 
(see Gibb et al. 2004 for a review). Hence, CO2 should be initially present in the vapor 
phase of the solar nebula, with C0 2 :CO = 1 to 4, a range of values that covers the ISM 
measurements (Gibb et al. 2004). 

Following the conclusions of Mousis et al. (2002), we also adopt for N2:NH3 a ratio 
close to the ISM value. However, even though current chemical models of the solar nebula 
predict that N 2 is more abundant than NH 3 , no observations have been published that would 
constrain this ratio further. In fact, some local processes may accelerate the conversion of 
N 2 to NH 3 in the solar nebula. This is for example the case if there exist some Fe grains that 
may have a catalytic effect on this conversion (Fegley 2000). Given these uncertainties, we 
compute the expected abundances of volatiles trapped in icy planetesimals using N 2 :NH 3 = 
1 and 10. 

Finally, it is important to note that since the evolution of the nebula proceeds from high 
to low temperatures, C0 2 crystalizes as a pure condensate prior to be trapped by water to 
form a clathrate hydrate. Therefore, we assume that solid C0 2 is the only existing condensed 
form of C0 2 present in the solar nebula. 

With these assumptions, Eq. 1 and 2 together with the thermodynamical evolution of 
the disk from A05 allow us to calculate the mass ratio of trapped volatiles to water, in the 
icy planetesimals accreted by the growing Jupiter (see Table 2). 

3.2. Enrichments in volatiles in Jupiter 

Using the results presented in Table 2, we can calculate the enrichments in volatiles 
in Jupiter. Since the composition of planetesimals does not vary significantly along the 
migration path of Jupiter (between 15 AU to 5 AU, see Table 2), we assume that all the ices 
of accreted planetesimals have an identical composition. Furthermore, in our definition of 
the ices/rocks ratio (I/R), we include all the types of ices (H 2 0, C0 2 , hydrates and clathrate 
hydrates) . 

For the sake of comparison, we adopt first the same initial ratios for C and N species in 
vapor phase as in Hersant et al. (2004), namely C0 2 :CO:CH 4 = 0:10:1, and N 2 :NH 3 = 10. 
With these values, we reproduce their fit to the observed enrichments in volatiles. However, 
in order to obtain these values, 36.3M® of water and 9.7M e of other volatiles have to be 



- 6- 



accreted. Regardless of the assumed value for the planetesimals I/R ratio, this yields a total 
amount of heavy elements significantly exceeding the upper limit derived by G04. We note, 
as did Hersant et al. (2004), that the calculated enrichment of sulfur is significantly higher 
than the observed value. This problem stems from the assumption that S was only present 
under the form of H2S in the solar nebula. However, owing to the production of sulfur 
compounds from H2S at high temperature and pressure early in the nebula, the abundance 
of H 2 S may be subsolar by the time planetesimals form (Hersant et al. 2004). 

If we now assume C02:CO:CH4 = 30:10:1 and N2:NH3 = 1 (our nominal model) in 
vapor phase in the solar nebula, we also reproduce the enrichments in volatiles (see Table 

3) . However, the minimum mass of accreted ices is found to be 25.4M e , including at least 
11.3M e of water and 14.1M e of other volatiles. The resulting total content in heavy elements 
in Jupiter is therefore compatible with G04, provided the I/R ratio of accreted planetesimals 
is greater than ~ 1.5. This also translates into a minimum O/H value in Jupiter's envelope 
equal to ~ 4.2 x 10~ 3 (ie ~ 5 times the solar value of AG89) or about half the value obtained 
by Hersant et al. (2004). Finally, we note that the amount of water needed corresponds, 
if vaporized, to an abundance relative to H2 of at least ~ 9.5 x 10~ 4 , or ~ 2 times its 
gas-phase abundance in the solar nebula (see Table 1). Again, we explain this difference by 
sedimentation and drift effects of icy grains (Supulver and Lin 2000). 

We have estimated the sensitivity of our results to changes in the initial C02:CO:CH 4 
and N 2 :NH 3 ratios. Keeping N 2 :NH 3 = 1, we find that it is possible to fit the measured 
enrichments while remaining compatible with G04 for C02:CO = 1,2,3 or 4 (see Table 4). 
However, if N 2 :NH 3 is increased to 10, only C02:CO=l gives surface abundances compatible 
with the measured enrichments and a total amount of heavy elements compatible with G04 
but only provided the mean I/R ratio of planetesimals is at least equal to ~ 4.4 (see Table 

4) . 

Recently, the solar abundances have been subject to some revisions (Allende Prieto et 
al. 2001, 2002; see Lodders 2003 (L03 in the following) for a review). With these new values, 
the enrichments of C, N, S, Ar, Kr, and Xe are now 3.7 ± 1.0, 5.5 ± 2.4, 2.2 ± 0.3, 2.2 ± 0.8, 
2.0 ± 0.7, and 2.0 ± 0.5, taking into account the uncertainties on the solar abundances of 
L03. These enrichments can still be accounted for in our nominal model using the new solar 
abundances, except for N which is found to be enriched only by a factor 2.5 over the value 
given gy L03. However, note that N/H in Jupiter has been recently revised down (Wong et 
al. 2004), the resulting enrichment of N (relative to the values of L03) being 3.9 ± 1.9. In 
that case, our nominal model is compatible with all the abundances of the afore-mentionned 
elements in Jupiter. Finally, this requires the accretion of 20.6M© of ices, and translates 
into a minimum O/H in Jupiter's envelope equal to ~ 3.3 x 10~ 3 . In order to be compatible 
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with G04, the I/R ratio of accreted planetesimal must be greater than ~ 1. 

4. Summary and discussion 

In this Letter, we have presented calculations of the enrichments in volatile species in 
Jupiter's atmosphere using the latest giant planet formation models of A05 which include 
migration and disk evolution. The computed enrichments are therefore fully consistent with 
the thermodynamical evolution of the disk and with the growth of the giant planet. 

In our nominal model (C0 2 :CO:CH 4 = 30:10:1, and N 2 :NH 3 =1), the enrichments in 
volatiles observed in Jupiter's atmosphere result from the accretion of 25.4 M e of icy plan- 
etesimals. Assuming that these planetesimals contain at least ~ 60% of ices (I/R > 1.5), 
the total amount of heavy elements accreted is compatible with the determinations of G04. 
This consistency between surface abundance measurements and total content in heavy ele- 
ments remains valid for a range of C0 2 :CO and N 2 :NH 3 ratios albeit requiring in some cases 
stronger constraints on the minimum I/R ratio of accreted planetesimals (from ~ 1.4 to 
~ 4.4). We note that if we ignore the presence of pure C0 2 ice but assume that C is trapped 
in form of CO and CH4 only, we were capable to match the measured enrichments while 
remaining compatible with G04 only if we assumed the accretion of pure ice planetesimals 
(I/R = 00). 

Furthermore, our nominal model predicts that the O/H ratio in the atmosphere of 
Jupiter must be at least ~ 4.2 x 10~ 3 which is half the one obtained by Hersant et al. 
(2004). This ratio can possibly be measured by future space missions. Interestingly, the 
minimum I/R ratio (~ 1.5) required to satisfy the upper limit derived by G04 happens to be 
consistent with the estimated ices content of the Galilean icy satellites (50 to 60%, see Sohl 
et al. 2002). Our results are then compatible with the formation of regular icy satellites from 
the accretion of planetesimals originating from the nebula without having been vaporized 
inside the subdisk (Mousis & Gautier 2004, Mousis et al. 2005). 

We conclude that the enrichments in volatiles in Jupiter's atmosphere can then be 
explained in a way compatible with both the internal structure models of the giant planet, 
and the average composition of its regular icy satellites. Furthermore, if one assumes a value 
of I/R in planetesimals of the order of ~ 1.5, the enrichments in volatiles are also compatible 
with the high mass model of G04. This model is however characterized by a low core mass 
while the model of A05 predicts a core mass of ~ 6M ffi . Thus, we must assume that the core 
has been eroded over time (see G04 for details) and that the volatiles which have not been 
released during the accretion phase itself have been released in the atmosphere during this 
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process. 
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Table 1: Gas-phase abundances (molar mixing ratio with respect to H 2 ) of major species in 
the solar nebula (from AG89) for C0 2 :CO:CH 4 = 30:10:1 and N 2 :NH 3 = 1 (our nominal 
model). 



Species i 






Species i 









1.71 x 10" 


-3 


N 2 


7.47 x 10- 


-5 


c 


7.26 x 10" 


-4 


NH 3 


7.47 x 10- 


-5 


N 


2.24 x 10- 


-4 


S 


3.24 x 10" 


-5 


H 2 


4.86 x 10" 


-4 


Ar 


7.26 x 10" 


-6 


co 2 


5.31 x 10" 


-4 


Kr 


3.39 x 10" 


-9 


CO 


1.77 x 10" 


-4 


Xe 


3.39 x 10- 


10 


CH 4 


1.77 x 10- 


-5 
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5 AU 15 AU 



C0 2 :H 2 


9.43 x 


10" 


-1 


9.28 x 


10" 


-1 


CO:H 2 


1.47 x 


10" 


-1 


1.48 x 


10" 


-1 


CH 4 :H 2 


9.44 x 


10' 


-3 


9.39 x 


10" 


-3 


N 2 :H 2 


6.07 x 


10" 


-2 


6.11 x 


10" 


-2 


NH 3 :H 2 


5.45 x 


10" 


-2 


5.34 x 


10" 


-2 


H 2 S:H 2 


3.16 x 


10" 


-2 


3.11 x 


io- 


-2 


Ar:H 2 


6.83 x 


10" 


-3 


6.95 x 


10" 


-3 


Kr:H 2 


8.01 x 


10" 


-6 


8.10 x 


10" 


-6 


Xe:H 2 


1.52 x 


10" 


-6 


1.51 x 


10" 


-6 



Table 2: Calculations of the ratios of trapped masses of volatiles to the mass of H 2 ice (Yi in 
the text) in planetesimals formed at 5 and 15 AU in the solar nebula. Gas-phase abundance 
of H 2 is equal to 9.5 x 10~ 4 (~ 2 times the value quoted in Table 1, see text), and gas-phase 
abundances of elements are assumed to be solar (AG89) with C0 2 :CO:CH 4 = 30:10:1 and 
with N 2 :NH 3 = 1 in vapor phase in the solar nebula. The abundance of H 2 S is subsolar (see 
text). 
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Species 


Observed 


Calculated 


Ar 


2.5 ±0.5 a 


2.15 


Kr 


2.7 ±0.5 a 


2.30 


Xe 


2.6 ±0.5 a 


2.80 


C 


2.9 ±0.5 b 


3.10 


N 


3.6 ±0.8 b 


2.80 


S 


2.5 ±0.15 b 


2.35 (3.40 without correction) 



Table 3: Observed enrichments in volatiles in Jupiter, and calculated enrichments in volatiles 
in the nominal model (C0 2 :CO:CH 4 =30:10:1 and N 2 :NH 3 = 1). The observed values are 
taken from Mahaffy et al. (2000) (a) and Atreya et al. (1999) (b), and are relative to solar 
abundances of AG89. The calculated abundance of sulfur (3.4 without correction) is revised 
down to fit the observed value (H 2 S/H 2 = 0.7 x (S/H 2 ) ; see text). 
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N 2 :Nfh 


t = 1 




N 2 :NH 3 = 10 




C0 2 :CO = 4 


C0 2 :CO = 3 


C0 2 :CO = 2 


C0 2 :CO = 1 


C0 2 :CO = 1 


M water /M e 


9.9 


11.3 


13.5 


17.9 


20.8 


M ices /M e 


24.4 


25.4 


26.9 


30.0 


34.2 


I/R 


1.4 


1.5 


1.8 


2.5 


4.4 


0/H 


4.7 


5 


5.3 


6.1 


7 



Table 4: Minimum mass of accreted water (M water ), minimum mass of accreted ices (Mi ces ), 
and minimum I/R ratio required to fit both the observed enrichments in volatiles, and the 
constraints of G04. The last line gives the minimum O/H abundance (compared to solar 
value of AG89) in the atmosphere of Jupiter. 
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Log 10 [Total pressure (bars)] 



Fig. 1. — Stability curves of the condensates considered in the present work, and evolution- 
ary tracks of the nebula at 5 and 15 AU. Abundances of various elements are solar, with 
C02:CO:CH 4 = 30:10:1 and N 2 :NH 3 = 1 in vapor phase. The condensation curve of CO2 
pure condensate (solid line) is shown together with that of the corresponding clathrate hy- 
drate (dashed line). Species remain in the vapor phase as long as they stay in the domains 
located above the curves of stability. 



